17 Intermittency (and quasiperiodicity) 

In this lecture we discuss the other two generic routes to chaos, intermittency 
and quasiperiodicity. 

Almost all our remarks will be on intermittency; we close with a brief de- 
scription of quasiperiodicity. 

Definition: Intermittency is the occurrence of a signal that alternates ran- 
domly between regular (laminar) phases and relatively short irregular bursts. 

In the exercises we have already seen examples, particulary in the Lorenz 
model (where it was discovered, by Manneville and Pomeau). 

Examples: 

• The Lorenz model, near r = 166. 

Figure la,b Manneville and Pomeau (1980) 

• Rayleigh-Benard convection. 
BPV, Figure IX.9 



17.1 General characteristics of intermittency 

Let r = control parameter. The following summarizes the behavior with 
respect to r: 

• For r < rj, system displays stable oscillations (e.g., a limit cycle). 

• For r > ri (r — ri small), system in in the intermittent regime: stable 
oscillations are interrupted by fluctuations. 

• As r — > T{ from above, the fluctuations become increasingly rare, and 
disappear for r < T{. 

• Only the average intermission time between fluctuations varies, not their 
amplitude nor their duration. 
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We seek theories for 



• Linear stability of the limit cycle and "relaminarization." (i.e. return to 
stability after irregular bursts). 

• Scaling law for intermission times. 

• Scaling law for Lyaponov exponents. 



17.2 One-dimensional map 

We consider the instability of a Poincare map due to the crossing of the unit 
circle at (+1) by an eigenvalue of the Floquet matrix. 

This corresponds to the specific case of Type I intermittency. 

Let u be the coordinate in the plane of the Poincare section that points in 
the direction of the eigenvector whose eigenvalue A crosses +1. 

The lowest-order approximation of the 1-D map constructed along this line 
is 

u' = \{r)u. (39) 
Taking A(r^) = 1 at the intermittency threshold, we have 

u' = X(n)u = u. (40) 

We consider this to be the leading term of a Taylor series expansion of u'(u, r) 
in the neighborhood of u = and r = r^. 

Expand to first order in (r — r^) and second order in u: 

dv! 



dv! 

u(u, r) ~ w'(0, Tj) + u 



du 



1 o 6>V 



+ (r - Ti) 



0,r l 



Evaluating equation (39), we find that the first term vanishes: 

u'(u = 0, r = Vi) = 0. 
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From equation (40), we have 



dvf_ 
du 



= Kn) = 



0,n 



Finally, rescale u such that 



1 6>V 



2 du 2 



= 1 



and set 

The model now reads 



e oc {r — Vi). 



u' = U + £ + W 2 . 

where £ is now the control parameter. 



Graphically, we have the following system: 



• £ < 0, i.e. r < T{. 

• w_ is stable fixed point. 

• u + is unstable. 



• £ = 0, i.e. r = r-j. 

• u' is tangent to indentity map. 

• w_ = m + = is marginally stable. 
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• e > 0, i.e. r > Vi. 

• no fixed points. 



For e < 0, the iterations look like 




• U- is an attractor for initial conditions 
u < u + . 

• For initial conditions u > u+, the itera- 
tions diverge. 
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The situation changes for £ > 0, i.e. r > T{ 



• No fixed points. 

• Iterations beginning at u < drift towards 
u > 0. 




The fixed points of u'(u) represent stable oscillations of the continuous flow. 

Thus for u ~ 0, the drift for £ > corresponds to a flow qualitatively similar 
to the stable oscillations near u = for £ < 0. 
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However, when e > 0, there is no fixed point, and thus no periodic solution. 

The iterations eventually run away and become unstable — this is the inter- 
mittent burst of noise. 

How does the laminar phase begin again, or "relaminarize" ? 
Qualitatively, the picture can look like 



Note that the precise timing of the turbulent burst is unpredictable. 

The discontinuity is not inconsistent with the presumed continuity of the 
underlying equations of motion — this is a map, not a flow. 

Moreover the Lorenz map itself contains a discontinuity, corresponding to the 
location of the unstable fixed point. 

17.3 Average duration of laminar phase 

What can we say about the average duration of the laminar phases? 
Writing our theoretical model as a map indexed by k, we have 




u 



U k +1 = U k + £ + u\. 
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For Uk+i — Uk, we can instead write the differential equation 

du 9 

The general solution of this o.d.e. is 

u(k) = e 1/2 tan \e 1/2 (k - k ) 

Take ko = 0, the step at which iterations traverse the narrowest part of the 
channel. 

We thus have 

u{k) = e 1/2 taa (V /2 fc). 
We see that u(k) diverges when 

£ V*k = ±- or k = ±-e- 1 ' 2 . 
2 2 

The divergence signifies a turbulent burst. 

When k ~ £~ 1//2 , itfc + i — is no longer small, and the differential approxi- 
mation of the difference equation is no longer valid. 

Thus: if r = time (oc number of iterations) needed to traverse the channel, 
then 

t oc e~ 1/2 or t (X (r - r z )~ 1/2 . (41) 

Thus the laminar phase lasts increasingly long as the threshold r = T{ is 
approached from above. 



17.4 Lyaponov number 

We can also predict a scaling law for the Lyaponov number. 

Near the fixed point (u ~ 0, e > 0), the increment Siik+i due to an increment 
Uk is, to first order, 

5u k+ i ~ XiSu k 
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where Ai is eigenvalue that passes through (+1). 
After N iterations, 

5u N ~ AivAiv-iAiv_2 • • • \\8u\. 
Suppose N ~ the duration of the laminar phase. Then 

Xn > 1 and \n-i — A^v-2 — • • • — Xi — 1. 

The Lyapunov number A is 

i 

where the last relation USed equation (41). (Recall that In A = Lyaponov exponent.) 

Results from the Lorenz model verify this prediction. The "intermittent 
channel" of the Lorenz map is seen in 

BPV, Figure IX. 14 

and the associated £ 1//2 scaling of the Lyaponov number is seen in 
BPV, Figure IX. 15 

Behavior qualitatively similar to that predicted by our model has been ob- 
served in the B-Z reaction: 

BPV, Figure IX. 16-17 
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17.5 Quasiperiodicity 

Finally, we make a few remarks about the third universal route to chaos, 
known as quasiperiodicity. 

Recall that there are 3 generic ways in which a limit cycle on a Poincare map 
may become unstable: An eigenvalue A of the Floquet matrix (the Jacobian 
of the map) crosses the unit circle at 

• +1 (as in the example of intermittency above); 

• —1 (as we saw in the introduction to period doubling); and 

• A = a ± i/3, |A| > 1. This corresponds to the transition via quasiperiod- 
icity. 

As we have seen, the latter case results in the addition of a second oscillation. 

This is a Hopf bifurcation: the transformation of a limit cycle to a quasiperi- 
odic flow, or a torus T 2 . 

The route to chaos via quasiperiodicity describes how a torus T 2 (i.e., a 
quasiperiodic flow) can become a strange attractor. 

17.5.1 An historical note 

In 1944, the Russian physicist Landau proposed a theory for the transition 
from laminar flow to turbulence as the Reynolds number is increased. 

Briefly, he envisioned the following sequence of events as Re increases beyond 
Re c : 

• Laminar flow (constant velocity) becomes periodic with frequency fi by 
a Hopf bifurcation. 

• Period flow^quasiperiodic flow; i.e., another Hopf bifurcation. The sec- 
ond frequency is incommensurate with j\. 
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• More incommensurate frequencies f%, f<±, . . . , f r appear in succession (due 
to more Hopf bifurcations). 

• For r large, the spectrum appears continuous and the flow (on a torus 
T r ) is aperiodic (i.e., turbulent). 

Recall that we have learned previously that, for dissipative flows, 

dimension of phase space > attractor dimension. 

Thus a consequence of Landau's theory is that a system must have many 
degrees of freedom to become chaotic. 

We now know, however, from the work of Lorenz, that 

• 3 degrees of freedom suffice to give rise to a chaotic flow; and 

• the chaos occurs on a strange attractor, which is distinct from a torus 
(since trajectories diverge on the strange attractor). 

17.5.2 Ruelle-Takens theory 

Lorenz 's observations were deduced theoretically by Ruelle and Takens in 
1971. 

The Ruelle-Takens theory is the quasiperiodic route to chaos. As a control 
parameter is varied, the following sequence of events can occur: 

• Laminar flow —> oscillation with frequency f\. 

• A second Hopf bifurcation adds a second (incommensurate) frequency 

• A third Hopf bifurcation adds a third frequency fa. 

• The torus T 3 can become unstable and be replaced by a strange attrac- 
tor. 
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The transition is demonstrated beautifully in terms of changing power spectra 
in the Rayleigh-Benard experiment described by 

Libchaber et al., Figure 15 

Libchaber, A., Fauve, S. and C. Laroche.1983. Two-parameter study of the routes to chaos. Physica 
D. 7: 73-84. 

Note that the Rayleigh number of the two spectra varies by less than 1%. 

Such a transition can also be seen in Poincare sections, such as the Rayleigh- 
Benard experiment of 

BPV, Figures VII.20, VII.21 
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